Domain wall dominated universes 
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We consider a cosmogony with a dark matter component consisting of a network of frustrated 
domain walls. Such a network provides a solid dark matter component with p — — (2/3)p that 
remains unclustered on small scales and with f2dw ~ 0.7 can reconcile a spatially flat universe with 
the many observations indicating i}^ « 0.3. Because of its large negative pressure, this compo- 
nent can explain the recent observations indicating an accelerating universe without recourse to a 
non- vanishing cosmological constant. We explore the viability of this proposal and prospects for 
distinguishing it from other kinds of proposed dark matter with significant negative pressure. 



Recent observations of apparent luminosities of Type la supernovae (SNIa) at moderate redshift ((z) « 0.6) suggest 
that the expansion of the universe is now accelerating, indicating a form of matter with a significant negative pres- 
'. sure From a theoretical standpoint, these observations are remarkable because most conceivable contributions to 
' the stress-energy entail a positive (radiation) or vanishing (matter) pressure, rather than a large negative pressure. 
Two proposals to explain these observations are a non-vanishing cosmological constant (A) or a very slowly rolling 
scalar field, often dubbed quintessence. Both proposals, however, are plagued with formidable fine tuning problems, 
psj The fine tuning problems of A 7^ are well known Q . For the alternative of a slowly rolling scalar field |^ extraordi- 
^ narily flat potentials are required, so that the field is unable to roll to the true minimum by the present day. If the 
. necessary flatness of such a potential is characterized by a mass scale to, one requires to ~ 10"^"^ eV. 
I In this letter we suggest another form of dark matter with significant negative pressure — a solid dark matter (SDM) 
QQ - component with the properties of a relativistic solid The term solid here denotes a substance with a harmonic, 
non-dissipative resistance to pure shear (volume preserving) deformations. A perfect fluid with negative pressure 
0^ {w = p/p < 0) would have an imaginary sound speed, indicating instabilities most severe on the smallest scales. 
For a solid, however, with a sufficiently large shear modulus these instabilities are removed. Since the sound speed 
. of perturbations of the solid should comprise a substantial fraction of the speed of light, its Jeans length today is 
Oh' comparable to the size of the current horizon, and hence an SDM component would remain unclustered on the smaller 
Q . scales over which fl^ is measured and thus evade detection. 

i—^ ' To compute the effect of an exotic dark matter component on the evolution of cosmological perturbations, more is 
c/3 . required than just knowing how w evolves throughout cosmic history ||^. As we will see, it is possible to construct 

■ dark matter components for which w agrees at all times, but which respond differently to cosmological perturbations. 
^ ,. Because of the non-vanishing shear modulus in the SDM case, large anisotropic stresses are generated, whereas in 

models with a A 7^ or a slowly rolling scalar field these stresses vanish. Furthermore, long-wavelength gravity waves 
entering the horizon at late times acquire an effective mass because of the energetic cost associated with pure shear 

■ deformations of spacetime. In a previous paper Q two of us (MB and DNS) developed a formalism for computing the 
evolution of cosmological perturbations in the presence of a SDM component and computed the cosmic microwave 
background (CMB) anisotropy on large scales. Here we discuss the viability of such scenarios for creating cosmic 
structure and the prospects for distinguishing SDM from other types of dark matter with significant negative pressure 
using future measurements of the CMB. A subsequent paper |^ will elucidate in more detail the differences between 
these scenarios and other models, such as quintessence and A 7^ 0. 

A possible microphysical origin for SDM is a frustrated network of topological defects, either cosmic strings or 
domain walls (see ref. ^ for a review of cosmic defects). The standard picture for the formation and evolution of 
topological defects is that of a random network of defects produced in a cosmological phase transition which then 
evolves toward a self-similar scaling regime with the number of defects within a Hubble volume approaching a fixed 
number However, such behavior relies on the ability of the defects to untangle and lose energy essentially as 
fast as allowed by causality, and in more complicated models this need not occur. In particular, in theories with 
several species of defects or with non-Abelian symmetries, topological obstructions to such untangling may arise 
and defect-dominated evolution can occur. The basic features of such evolution would be that after the decay of an 
initial transient the number of defects per co-moving volume approaches a constant ||l^,0. For the simple domain 
walls of most interest here this implies an equation of state with w = —2/3 (pdw ^ o^^) and for strings w = —1/3 
(Pstr ~ a^^)) although other values of w may be possible as well. Simulations of non-Abelian cosmic strings |ll| ] 
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FIG. 1. Rescaled CMB Moments for SDM Models. In the first four panels £{£ + l)Ci versus £ is plotted for 
w = —1/3, —1/2, —2/3, —5/6 for several values of the longitudinal sound speed for a universe with Qrn ~ 0.3, fiadm = 0.7. 
The solid curve corresponds to Cs ~ 0.0 and Cs = 0.2, 0.4, 0.6, 0.8 are indicated by the line styles dotted, dashed, dot-dash, 
and dot-dot-dash, respectively, rescaled so that the horizon size at recombination subtends the same angle in the sky today 
as for a. w — ~1 model with the same f2m. The bottom middle panel indicates the rescaled and renormalized moments for the 
quintessence models with w — —1/3, —1/2, —2/3, —5/6, and the bottom right panel shows all models of the previous panels for 
comparison. 



indicate behavior leading to a string-dominated universe and simulations of domain walls find that the density of 
walls in a model with several vacua falls less rapidly than would be required by scaling. In axion models, where an 
additional symmetry breaking U{1) Zjq allows for N domain walls to meet at a string, scaling behavior has been 
observed p^ . The precise criterion for when defect domination takes place is not clear at present, although it seems 
clear that such models can exist. This is currently under investigation. 

The symmetry breaking scale r\ required for the formation of a network of frustrated domain walls with f2dw ~ 1 
today is around r] ~ 100 KeV, assuming a phase transition at T « 77, with initially one wall per horizon volume 
and a network immediately settling down to an equilibrium configuration which is subsequently swept along by the 
Hubble flow. This estimate, however, is subject to considerable uncertainties in either direction. For example, a 
long transient before settling down to an equilibrium conflguration could considerably raise ry. The mean separation 
between walls is approximately 30 parsecs, much smaller than the scale on which the network clusters in response 
to gravitational perturbations. It is therefore justified to treat the network as a continuum solid when studying its 
response to cosmological perturbations. 

In a recent article it was suggested that a quintessence model with a flat geometry w w —2/3, ilc = 0.25, 
fib = 0.05, h = 0.65, n = 1 gave a adequate flt to data probing a range of epochs and scales, where fix is the 
fractional density relative to the critical density in particle species X (X=c is CDM, X=b is baryons), the Hubble 
constant is Ha = lOO/i km see" ^ Mpc~^ and n is the spectral index of the initial density perturbations. One might, 
therefore, wonder whether a domain wall dominated model with fidw ~ 0.7 might fit the data equally well and, in 
fact, a number of the calculations presented in ref. apply equally well to the case of a SDM model. However, 
important differences in how perturbations evolve once the SDM component comes to dominate must be considered 
for an accurate comparison to the data. Due to the apparently good fit of the w = —2/3 case, most of the discussion 
focuses on the domain wall dominated case, but we also discuss the string dominated case and other values of w. 
Except where expressly stated, all our results use the parameter choices listed above. 

We have included the evolution of a SDM component and its effects on the other perturbation variables into the 
standard Einstein-Boltzmann solver CMBFAST. As well as specifying w, this requires the introduction of a another 
parameter Cg which is the sound speed of scalar perturbations in the solid, related to the the vector sound speed, Cu, 
by Cj = 4c^/3 -I- w. For a given w, the evolution of the Newtonian potentials $ and ^' at late times differs as one 
varies c^, and also differs from quintessence models. This leads to distinct integrated Sachs- Wolfe (ISW) contributions 
to the CMB anisotropics. At recombination, when the density of the SDM component is negligible (for w <, —1/3), 
these differences are suppressed and the small-angle anisotropics are identical for the same initial fluctuations. For 
this reason, when the results of our computations are plotted in Fig. they are normalized at an angular scale 
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TABLE I. The distinguishability of various types of solid dark matter and quintessence using the CMB. The table above 
indicates the expectation value of the natural logarithm of the likelihood of the model corresponding to the column relative to 
that of the row assuming that the CMB moments were generated by the model corresponding to the column. For all models 
Qm ~ 0.3. As described in the text, the CMB multipole moments of the models were rescaled so that the overall amplitudes and 
the shapes coincide at large £. These relative likelihoods take into account only cosmic variance and assume full sky coverage. 
Observational errors and the fact that varying other cosmological parameters affects the shape of the low-^ moments is not 
taken into account. 



corresponding to £ = 500 and spectra for different values of Cs are then identical for i > 100. Models with different 
w, however, have different peak positions with all the other cosmological parameters fixed, since the angular diameter 
distance £d — fcp'7o depends on w. The wavenumber kp corresponds to the first acoustic peak [ p5[ and any change in 
770 due to a variation in w can be offset by modifying kp, as discussed in ref. ||l6|. In Fig. |^, the angular scale has 
been rescaled so that all the models have the angular diameter distance of the model with w = ~1. Distinguishing 
among the range of models considered will require accurate measurements of the CMB anisotropics on large scales. 

Table 1 examines the significance of these differences by comparing flat models with identical cosmological param- 
eters, but with differing properties for the component with significant negative pressure rescaled and normalized as 
described above. An estimate of the significance may be obtained by ignoring observational uncertainties and incom- 
plete sky coverage because of the galaxy and considering only the effect of cosmic variance, arising from the fact that 
for each £ one is able to observe only (2£+ 1) realizations of a Gaussian random process. Assuming that a model A is 
correct, the expectation value of the natural logarithm of the relative likelihood of model A relative to a model B is 
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For a domain wall dominated universe [w — —2/3) we find it virtually impossible to distinguish between differing 
sound speeds or between SDM and quintessence, while for w = —1/3 we observe substantial differences between 
differing sound speeds, and between SDM and quintessence. Moreover, all the models are distinguishable from A 
model with the same value of il,„. The relative likelihood of different models increases as Vtya is decreased 

An important test for any cosmological model is whether it can create the large-scale structure (LSS) seen today 
when normalized to the the CMB anisotropics observed by COBE. The simplest version of this is to compare the 
COBE normalized value of the variance of the density field in spheres of radius 8/i~"'^Mpc, with that obtained 
from the observations of X-ray clusters ||l^. Models with w = —1 pass this test, but two effects can cause models 
with w > — 1 to predict a lower erg when normalized to COBE: a large integrated Sachs- Wolfe component reduces 
the predicted value of cts, and for larger w the growth of perturbations becomes stunted earlier because the SDM 
(or quintessence) component begins to dominate the universe earlier. To investigate the viability of these models we 
compute the COBE normalized erg for w = —2/3 and a range of values in the — h plane, which are then compared 
to the recently computed observational value crg^^^i = (0.5 -0.16) ±0.1, where 7 = 0.21 -0.22w + 0.33f^ni + 0.259 
and Q = {n — 1) + [h — 0.65), valid in the range of parameters considered. Fig. || shows the regions of fim — h plane 
which pass this test for spectral indices n = 1.0 and n = 1.1. For n = 1.0, the preferred values for quintessence models 
h = 0.65 and — 0.3 are marginally incompatible with the quoted values for (Tg, but this situation is easily rectified 
by increasing h. Furthermore, these parameters are compatible with n = 1.1, although the viability of scenarios with 
w — —1/3 would require even larger values of n. The constraint from the upper bound on erg is weaker because of the 
possibility of including a tensor contribution to the CMB anisotropics . 
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FIG. 2. Domain wall dominated models compatible with constraints on ag from X-ray clusters. The solid line indicates the 
stringent lower limit and the dotted line the upper limit, which can be weakened by including a tensor component. The vertical 
and horizontal axes are h and $lm, respectively. The left and right panels indicate the spectral indices n = 1.0 and n — 1.1, 
respectively. In both cases the longitudinal sound speed is Cs = 0.2. 



We have shown that the SDM proposal, and in particular a domain wall dominated universe, is compatible with LSS 
via the COBE normalization/as test and that the inclusion of a SDM component, as opposed to quintessence, can have 
some potentially interesting effects on the CMB anisotropics. If the SNIa results are confirmed, an important question 
will be how to distinguish between SDM and quintessence models. Some possibilities include cross-correlating the CMB 
with LSS using, for example the X-ray background |lj], or the gravitational lensing of the CMB jl^. Observations 
of gravitational lensing can measure the evolution of the power spectrum with redshift ||20| , |2l|] , and in particular the 
proposed Dark Matter Telescope would measure the growth of perturbations with high accuracy. 
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